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Abstract 


The  purpose  of  this  thesis  is  to  investigate  and 
analyze  bandpass  sampling  techniques.  Bandpass  sampling 
theory  is  shown  to  be  an  extension  of  the  Shannon  low-pass 
sampling  theory.  Five  bandpass  sampling  techniques  are 
investigated:  direct  sampling  (also  called  first-order 
sampling),  pth-order  sampling  (in  particular  second-order 
sampling  where  p=2),  conventional  quadrature  sampling, 
quadrature  sampling,  and  complex  sampling  (also  called 
Hilbert  transform  sampling). 

Direct  sampling  is  a  single  channel  approach  to 
bandpass  sampling  while  the  other  approaches  require  two 
channels.  The  two  channel  techniques  offer  a  sampling  rate 
reduction  of  up  to  one  half  over  direct  sampling,  but  there 
are  trade-offs  to  consider  when  choosing  a  two  channel 
technique.  It  is  shown  that  the  effects  of  random  phase 
error  can  render  the  conventional  quadrature  sampling 
approach  useless.  Timing  is  critical  with  quadrature 
sampling.  Reconstruction  error,  due  to  hardware  timing 
errors,  is  shown  to  be  directly  proportional  to  the  center 
frequency  and  is  a  function  of  the  bandwidth  of  the  signal. 
Plots  were  developed  to  show  the  effect  of  timing  on  a  given 
bandpass  signal  when  the  center  frequency  and  the  bandwidth 


are  known. 


(t,  I*.  gl  ’*1- 


I.  INTRODUCTION 


Problem 


The  purpose  of  this  thesis  is  to  investigate  and 
analyze  various  bandpass  sampling  techniques. 


Scope 

This  thesis  considers  direct  sampling  (also  called 
first-order  sampling),  2nd-order  sampling,  quadrature 
sampling,  and  complex  sampling  (sampling  using  Hilbert 
Transforms).  There  are  various  implementations  of 
each  technique  and  at  times  these  categories  overlap.  The 
details  will  be  presented  in  Chapter  II. 


Background 

In  a  digital  processing  receiver  the  received  signal  is 
sampled,  quantized,  and  processed  by  a  digital  processor. 

The  signal  may  be  sampled  at  any  point  in  the  RF  or  IF 
portions  of  the  receiver.  A  block  diagram  of  a  digital 
processing  receiver  is  shown  in  Figure  1. 


Figure  1.  Digital  Processing  Receiver 


For  this  thesis,  uniform  sampling  is  assumed  to  occur 
in  the  IF  stage.  The  resulting  sampled  spectrum  may  either 
be  a  real  bandpass  spectrum  or  a  complex  low-pass  spectrum 
depending  upon  the  sampling  technique  employed. 

If  a  band-limited  waveform  is  sampled  at  a  uniform  and 
sufficiently  high  rate,  the  original  waveform  can  be 
reconstructed  from  the  sampled  values  by  appropriate 
interpolating  functions.  Shannon  showed  that  the  minimum 
sampling  rate  is  equal  to  twice  the  highest  frequency 
component  of  the  waveform  (9:519).  However,  for  bandpass 
waveforms  like  the  one  shown  in  Figure  2,  it  is 
theoretically  possible  to  sample  at  a  rate  of  cr/n  samples/s 
or  twice  the  bandwidth  of  the  bandpass  waveform.  Thus, 
depending  upon  the  bandpass  sampling  technique  used,  the 
sampling  rate  for  a  bandpass  signal  can  be  much  lower  than 
twice  the  highest  frequency. 


Figure  2.  Typical  Bandpass  Waveform  (13:907) 


Direct  sampling,  as  shown  in  Figure  3,  involves 
multiplying  the  bandpass  waveform  with  a  train  of  uniformly 
spaced  pulses.  The  spacing  of  the  pulses  (or  the  sampling 
rate)  is  governed  by  the  bandwidth  of  the  bandpass  waveform. 
Special  consideration  needs  to  be  given  to  the  sampling  rate 
to  avoid  the  problem  of  aliasing.  It  will  be  shown  that  for 
direct  sampling  of  the  bandpass  waveform  the  minimum 
sampling  rate,  fs  samples/s,  ranges  from  twice  the  bandwidth 
to  four  times  the  bandwidth  and  is  a  function  of  the 
relationship  between  the  center  frequency  and  the  bandwidth. 
This  will  be  examined  further  in  Chapter  II. 
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Figure  3.  Direct  Sampling  (10:16) 

Pth-order  sampling  is  an  extension  of  direct  sampling. 
Where  direct  sampling  involves  a  single  train  of  uniformly 
separated  samples  of  the  bandpass  waveform,  pth-order 
sampling  involves  p  trains  of  samples  as  shown  in  Figure  4 
for  the  special  case  of  p=2.  The  trains  are  slightly 
displaced  in  time  but  occur  at  the  same  frequency  (10:40). 
The  sampling  rate  for  each  of  the  p  channels  is  fs/p 
samples/s  for  an  average  rate  of  fs  samples/s. 
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Figure  4.  Second-Order  Sampling  (13:905) 

Two  implementations  of  quadrature  sampling  will  be 
discussed.  The  first,  from  this  point  on  to  be  referred  to 
as  conventional  quadrature  sampling,  is  sketched  in  Figure  5 
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Figure  5.  Conventional  Quadrature  Sampling  (10:37) 
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Here,  the  bandpass  signal  is  preprocessed  by  employing 
quadrature  demodulation  prior  to  sampling.  In  one  channel 
x(t)  is  modulated  by  cosw0t  and  low-pass  filtered  to 
remove  the  sum  frequency  term  to  obtain  the  in-phase 
component,  p(t).  Similarly,  the  quadrature  component, 
q(t),  is  obtained  from  the  sine  channel.  Because  of  the 
preprocessing  of  x(t),  the  components  p(t)  and  q(t)  can  be 
sampled  according  to  the  low-pass  sampling  theorem.  The 
original  waveform,  x(t),  can  be  reconstructed  from  samples 
of  the  low-pass  components. 

The  second  Implementation  of  quadrature  sampling, 
from  here  on  referred  to  as  simply  quadrature  sampling,  wa 
introduced  by  O.D.  Grace  and  S.P.  Pitt,  and  is  an 
application  of  second-order  sampling  (pth-order  sampling 
where  p=2).  Quadrature  sampling,  pictured  In  Figure  6, 
requires  uniform  sampling  of  both  the  bandpass  signal  and 
its  quarter  wavelength  (based  on  nominal  frequency  wQ) 
translation,  each  at  a  common  sampling  rate  depending  upon 
the  exact  relationship  between  wD  and  a  (the  bandwidth). 


Figure  6.  Quadrature  Sampling 


When  the  intersample  spacing  is  properly  chosen,  the 
bandpass  signal  can  be  reconstructed  in  its  entirety  from 
knowledge  of  the  sample  values;  moreover,  with  quadrature 
sampling,  the  (low-pass)  in-phase  and  quadrature  components 
of  the  bandpass  signal  have  a  simple  explicit  representation 
in  terms  of  samples  of  the  original  bandpass  signal  without 
the  preprocessing  required  in  conventional  quadrature 
sampling.  It  will  be  shown  that  the  minimum  sampling  rate 
for  quadrature  sampling  can  be  as  low  as  o/n  (average) 
samples/s  (2:1659,1662). 

Another  possibility  for  maintaining  the  (average) 
sampling  rate  equal  to  its  minimum  value  a/u  samples/s  is 
to  sample  both  the  bandpass  waveform,  x(t),  and  its  Hilbert 
transform,  xH(t),  at  a  rate  of  a/2*  samples/s.  This 
approach,  pictured  in  Figure  7,  is  called  complex  sampling 
(10:37-38) . 
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Approach 

The  approach  to  this  thesis  was  basic:  analyze  the 
bandpass  sampling  techniques  from  a  theoretical  standpoint 
and  draw  conclusions  from  the  analysis. 

In  chapter  II,  a  theoretical  analysis  of  each  bandpass 
sampling  technique  will  be  presented. 

Results  of  the  analysis  will  be  presented  in  Chapter 

III . 

Finally,  In  chapter  IV,  conclusions  are  drawn  based  upon  the 
theoretical  analysis  presented  in  Chapter  II  and  the  results 
detailed  in  Chapter  III. 


II.  THEORETICAL  ANALYSIS 


Introduction 

The  conversion  of  signals  from  analog  to  digital 
representation  involves  two  operations:  sampling  and 
quantization.  With  regard  to  sampling,  one  of  the  principle 
requirements  is  that  it  be  sufficiently  fast  to  represent 
the  waveform  adequately  (4:91).  But  how  do  you  determine 
what  is  "sufficiently  fast?"  For  low-pass  signals  the 
question  is  easily  answered  with  the  Shannon  sampling 
theorem  (15:43).  For  bandpass  signals,  sampling  is  a  little 
more  complicated.  A  direct  application  of  the  Shannon 
theorem  will  work,  but  there  are  several  techniques 
available  that  allow  us  to  sample  bandpass  signals  at  rates 
much  lower  than  twice  the  highest  frequency  component  in  the 
signal.  The  purpose  of  this  thesis  is  to  investigate  and 
analyze  bandpass  sampling  theory  and  to  examine  the 
various  techniques  that  can  be  used  to  implement  bandpass 
sampling.  This  chapter  is  devoted  to  the  development 
of  the  bandpass  sampling  techniques  from  a  theoretical 
standpo i nt . 

The  best  approach  to  the  topic  of  bandpass  sampling  is 
to  begin  with  a  review  of  the  Shannon  sampling  theorem.  As 
mentioned,  this  theorem  primarily  concerns  low-pass  signals, 
but  its  results  easily  extend  to  the  bandpass  case.  In 
addition,  by  reviewing  low-pass  sampling  theory,  sampling 


issues  such  as  aliasing  and  fold-down  (to  be  discussed  later 
in  this  chapter  -  and  again  in  Chapter  III)  can  be  addressed 
early,  before  examining  the  individual  bandpass  sampling 
techniques.  It  will  be  shown  that  the  aliasing  and  fold¬ 
down  issues  are  basically  handled  the  same  whether  talking 
about  low-pass  sampling  or  bandpass  sampling. 

The  approach  in  this  chapter  will  be  to  consider  the 
low-pass  sampling  theorem  first.  That  result  will  be 
extended  to  the  bandpass  case.  Then  the  individual  bandpass 
sampling  techniques,  mentioned  briefly  in  chapter  one,  will 
be  discussed. 

Shannon's  Low-Fass  Sampling  Theorem.  Shannon's  theorem 
states  that  a  signal  whose  Fourier  Transform  is  zero  outside 
the  Interval  Jwj  >  B  as  in  Figure  8,  can  be  uniquely 
represented  by  a  set  of  samples  of  the  waveform  taken  at 
intervals  of  n/B  seconds;  i.e.,  the  original  waveform  can 
be  completely  reconstructed  from  its  samples.  The  sampling 
frequency,  fs  =  B/n  ,  is  twice  the  highest  frequency 
present  in  the  waveform  and  is  usually  referred  to  as  the 
Nyquist  frequency  (4:91). 


Since  this  sampling  theorem  is  so  significant  to  the 
understanding  of  bandpass  sampling  theory,  two  different 
proofs  will  be  presented  here  (15:43,46).  However,  the 
proofs  are  not  mathematically  rigorous,  and  for  our 
purpose,  only  the  key  concepts  as  they  apply  to  bandpass 
sampling  will  be  reviewed  here. 

Proof  1.  Assume  X(w)  is  only  non-zero  along  the 
finite  interval  -B  <  w  <  B  ,  and  expand  it  in  a  Fourier 
Series  as 

m 

X(w)  =  C  cnexp{jnt0w)  (1 

n=-« 

where , 

t„  =  2n/2B  (2 

fB 

cn  =  ( 1 / 2 B )  I  X( w)exp{ - jntDw)dw  (3 

J-B 

Also,  the  Fourier  inversion  formula  tells  us  that 

rB 

x(t)  =  (l/2it)l  X  ( w)exp{  jwt  )dw  (4 

J -B 

Comparing  Eq  (3)  with  Eq  (4)  we  find  that 

cn  =  (n/B)x(-nt0)  =  ( n/B ) x { -nu/B ]  (5 

This  says  that  the  cn  are  known  once  x(t)  is  known  at  the 
points  t  *  nw/B  .  Plugging  these  values  of  cn  into  the 
series  expansion  for  X(w),  Eq  (1),  we  find  that 


X  ( w) 


(6) 


=  ( n/B ) E  x(-nn/B)exp( Jwnn/BJ 

n  =  -® 


Eq  (6)  is  essentially  a  mathematical  statement  of  the 
sampling  theorem.  It  states  that  X(w)  is  completely  known 
and  determined  by  the  sample  values,  x(-mt/B).  These  are 
the  values  of  x(t)  at  equally  spaced  points  in  time.  We  can 
plug  this  X(w),  given  by  Eq  (6),  into  the  inverse  Fourier 
integral,  Eq  (4),  to  find  x(t)  in  terms  of  its  samples. 

This  will  yield  the  following  result: 

fB 

x ( -mt/B )exp{ jnwn/Biexp{ jwt )dw 
-B 

(7) 

or,  after  simplifying 


x(t)  *  E  x ( -nn/B ) [ sin ( Bt+nu ) / ( Bt+nn ) ]  (8) 

n  =  -® 

The  original  signal,  x(t),  may  be  exactly  reconstructed  by 
interpolating  samples  spaced  t  =  nn/B  apart.  In  other 
words,  if  we  have  a  signal  x(t)  with  no  spectral  energy 
beyond  limit  B,  all  information  about  the  signal  is 
contained  in  a  sequence  of  samples  taken  at  a  uniform  rate 
greater  than  B/n  samples/s.  Note  Eq  (8)  has 
the  form 


x(t)  «  E  x ( kTs ) h ( t-kT- ) 


(9) 


where  Ts  Is  the  period  between  samples  taken  periodically  at 


times  kTs  ,  x(kTs)  are  the  samples,  and  h(t)  is  the 
low-pass  filter,  or  the  sampling  (interpolation)  function 


h ( t )  =  (sinBt )/Bt 


Notice  that  h(t)  has  the  form  sinx/x  ,  whose  transform  is 
a  spectral  window  centered  about  zero.  (The  importance  of 
this  spectral  window  will  be  more  obvious  after  Proof  2.) 

Eq  (9)  will  be  referenced  throughout  this  chapter. 

Although  this  proof  was  straightforward  and  it  gives  us 
x(t)  in  terms  of  its  samples,  it  provides  little  insight 
into  the  sampling  process.  A  more  intuitive  approach  is 
considered  in  Proof  2,  which  follows. 

Proof  2.  Assume  once  again  that  x(t)  has  the 
transform  X(w)  which  is  band  limited  to  the  Interval 
-B  <  w  <  B  as  in  Figure  8.  Next  consider  the  product  of 
x(t)  with  a  periodic  train  of  delta  functions,  p(t),  where 


p  ( t )  =  E  d (t-nTs ) 
n  =  -• 


The  sifting  property  of  the  Impulse  function  tells  us  that 


E  x(t)d(t-nT_) 


E  x(nTs )« (t-nTs ) 


Using  this  property,  we  find  that  the  sampled  version  of 
x(t),  xs(t),  is  given  by 


•-'.v.v.v  vvy.\ivSv. 


Xs(t)  =  x(t)p(t) 


(13) 


xs(t)  =  £  x(t)«(t-nTs) 
n  =  -® 


(14) 


xs(t)  =  £  x  (nTs  )  <J  ( t-nTs  ) 
n  =  -» 


(15) 


Figure  9  illustrates  the  relationship  of  x(t)  and  xs(t); 
i.e.,  xs(t)  is  x ( t )  sampled  at  intervals  of  Ts. 

To  show  that  x(t)  can  be  recovered  from  xs(t),  we  can 
use  the  frequency  convolution  theorem  to  get 


xs(t)  =  x(t)p(t)  < - >  (l/2n)X(w)*P(w)  =  Xs(w)  (16) 


where 


*  denotes  convolution 


CO 

P(w)  =  ( 2n/Ts ) £  «(w-nws) 

n  =  -<D 


(17  ) 


and 


ws  =  2u/Ts  (18) 

Recall  that  convolution  with  a  delta  function  simply  shifts 
the  original  function.  Therefore, 

•  CD 

£  X(w)*fl(w-nws)  =  £  X(w-nws)  (19) 

n=-»  n=-® 


13 


Xs(w)  =  (l/2n)X(w)*(2n/Ts )C  fi(w-nws) 

n  =  -a> 


Xs(w)  =  (1/TS)£  X ( w-nws ) 


This  Xs(w)  is  sketched  in  Figure  9.  From  the  sketch  of 


Xs(w)  it  is  clear  that  X(w)  can  be  identified.  In 


practice,  we  need  only  filter  Xs(w)  with  the  ideal  low- 


pass  filter,  H(w),  shown  in  Figure  10,  to  recover  X(w). 


This  is  the  spectral  window  referred  to  in  Proof  1.  The 


window  function  passes  all  frequency  components  in  the 


interval  -B  <  w  <  B  and  rejects  all  others 


In  Figure  9f  the  assumption  was  made  that  the  "humps" 


in  the  figure  do  not  overlap.  The  sampling  rate  is  set  at 


twice  the  bandwidth,  or  fs  =  B/n  samples/s.  Figure  11  is  a 


sketch  of  Xs(w)  with  fs  >  o/n  samples/s  and  fs  <  cr/n 


samples/s.  To  prevent  overlapping  spectra  the  point  labeled 


"ws  -  B"  must  fall  to  the  right  of  that  labeled  "B."  If 


this  were  not  true  as  in  Figure  11c,  the  result  would  be 


undersampling,  or  aliasing.  When  aliasing  occurs. 


l.  s 


It  is  impossible  to  uniquely  recover  X(w)  by  ideal  filtering 
or  any  other  means.  This  restriction  can  be  rewritten  as 
follows 

vs-B  >  B  < - >  ws  >  2B 

Since 

vs  =  2n/Ts 

this  restriction  becomes 

Ts  <  n/B  (24) 

In  other  words  the  samples  cannot  be  spaced  any  further 
apart  than  ti/B  samples. 

Simply  stated,  the  Shannon  sampling  theorem  requires 
that  x(t)  be  sampled  at  a  rate  greater  than  twice  the 
highest  frequency  component  in  x(t)  in  order  to  be  able  to 
recover  the  original  signal  from  the  sampled  data  form. 

This  process  is  summarized  in  Figure  12. 

Aliasing.  Aliasing  plays  an  important  role  in  bandpass 
sampling  just  as  in  low-pass  sampling.  To  be  able  to 
recover  the  original  signal  from  the  sampled  data,  none  of 
the  shifted  spectral  components  can  overlap.  This  makes 
bandpass  sampling  theory  a  bit  more  complicated  because  now 
two  spectral  bands  are  involved  in  the  bandpass  case  (one  at 
+w0  and  one  at  -w0  as  in  Figure  2)  as  opposed  to  only  one  in 
the  low-pass  case  (from  -B  to  +B  as  in  Figure  8). 


(22) 


(23) 
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Aliasing  can  occur  if  either  the  signal  is  not 


bandlimited  or  if  the  sampling  rate  is  too  low.  Theoretic¬ 
ally,  if  a  signal  is  not  bandlimited,  there  is  no  way  to 
avoid  the  aliasing  problem  with  the  basic  scheme  just 
discussed.  However,  the  spectra  of  most  realizable  signals 
may  be  assumed  to  be  bandlimited.  Furthermore,  a  common 
practice  is  to  filter  the  continuous  time  signal  before 
sampling  to  ensure  that  it  does  meet  the  band-limited 
criterion.  These  anti-aliasing  filters  are  also  important 
in  preventing  the  possibility  of  noise  or  interference 
folding  into  the  spectrum  occupied  by  our  signal  (16:47). 

The  use  of  antialiasing  filters  is  best  illustrated 
with  an  example.  Assume  the  speech  signal  in  Figure  13  is 
confined  to  the  frequency  interval  between  0.3  and  3.0kHz. 
Additionally,  there  is  interference  at  6kHz.  Our  concern  is 
what  happens  to  the  6kHz  sine  wave  as  a  result  of  sampling? 
If  the  composite  signal  is  sampled  at  8kHz  (a  standard  in 
telecommunications ) ,  the  6kHz  noise  will  be  folded  to  (or 
replicated  about)  the  2kHz  position  as  in  Figure  14.  The 
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Figure  14.  Example  of  Noise  Folding  Down  (12:64) 

interfering  signal  cannot  be  removed  without  loss  of  a 
portion  of  the  original  speech  signal.  The  noise  could, 
however,  be  removed  prior  to  sampling  by  filtering  the 
speech  signal  (12:63-64).  This  antialiasing  process  serves 
two  important  functions:  it  artificially  bandlimits  the 
signal  being  sampled;  and,  it  filters  out  noise  which  could 
otherwise  fold  down  and  make  signal  recovery  impossible 
(14:52).  More  will  be  said  about  antialiasing  in  Chapter 
III. 

This  background  material  lays  the  nece33ary  foundation 
for  discussing  the  individual  bandpass  sampling  techniques. 
References  will  be  made  throughout  the  rest  of  this  chapter 
to  this  material:  in  particular,  Eq  (9)  will  be  referenced 
with  the  development  of  each  bandpass  sampling  technique. 

Bandpass  sampling  involves  the  uniform  sampling  of  the 
bandpass  signal  by  one  or  more  train  of  impulses.  In  each 
case,  equations  for  the  bandpass  signal  will  be  developed  in 
terms  of  its  uniform  samples.  The  minimum  possible  sampling 


[•] 


rates  and  the  conditions  for  the  minimum  sampling  rates 
will  also  be  developed. 

The  purpose  of  this  chapter  is  to  develop  these 
bandpass  sampling  techniques.  Techniques  are  compared  in 
Chapter  III. 

Direct  Sampling 

Intr oduct 1  on .  Direct  sampling  will  be  developed  from 
two  points  of  view:  the  conventional  approach,  referred  to 
simply  as  direct  sampling,  is  developed  using  an  intuitive 
approach  to  give  a  ''feel”  for  bandpass  sampling  theory;  the 
second  approach,  called  first-order  sampling,  is 
mathematically  developed  and  reduces  the  sampling  problem  to 
an  application  of  the  low-pass  sampling  theorem. 

Conventional  Direct  Sampling.  The  first  bandpass 
sampling  technique,  direct  sampling,  can  be  viewed  as  the 
result  of  the  product  of  a  signal,  x(t),  and  a  periodic 
pulse  train,  p(t),  just  as  in  the  case  of  low-pass  sampling. 
The  bandpass  signal  and  the  sampling  process  are 
demonstrated  in  Figure  15.  Recall,  when  sampling  bandpass 
waveforms  we  no  longer  have  just  one  spectrum  which  is 
periodically  replicated;  Instead,  referring  back  to  Figure 
15b,  there  is  a  positive  and  a  negative  spectrum  which  get 
replicated.  The  sampled  spectrum  of  the  negative  frequency 
portion  of  X(w),  X_(w),  is  sketched  in  Figure  15d,  and 
corresponds  to  a  periodic  shift  every  2nfs  rad/s. 

Similarly,  the  sampled  spectrum  of  the  positive  frequency 
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portion  of  X(v) ,  X+(w),  Is  sketched  In  Figure  15e.  Figure 

15f  shows  the  combined  result  of  the  positive  and  negative 
spectrums  being  replicated.  Reconstruction  of  X(w)  will  be 
obtained  by  ideal  bandpass  filtering  of  Xs(w)  over  the 
original  bandwidth  of  X(w)  as  long  as  the  sampling  rate  is 
such  that  none  of  the  spectral  replicas  overlap  in  the 
original  bandpass  region. 

Analysis.  A  closer  examination  of  the  sampled 
spectrum  provides  insight  in  determining  the  minimum 
sampling  rate  possible.  It  will  be  shown  that  the  Nyquist 
sampling  rate  for  a  bandpass  waveform  of  twice  the 
bandwidth,  or  a/n  samples/s,  can  be  realized  only  when  the 
signal  is  positioned  in  such  a  way  that  the  upper  cutoff 
frequency  (w0  +  a/2)  is  an  integer  multiple  of  a.  Otherwise 
the  required  sampling  rate  will  be  a  value  between  o/n  and 
2a/n  samples/s.  To  determine  the  minimum  sampling  rate  for 
these  two  cases  we'll  consider  each  separately. 

First,  consider  Case  1  in  Figure  16  where  wQ  +  a/2  =  kcr 
(k  =  integer).  Just  as  we  reconstructed  the  low-pass 
waveform  in  our  earlier  examples,  we  wish  to  be  able  to 
reconstruct  our  sampled  waveform  over  the  original  bandpass 
spectrum.  Notice  in  the  Case  1,  where  w0  +  a/2  =  kcr  and 
fs  =  a/n  samples/s,  none  of  the  replicated  spectrums 
overlap  the  original  bandpass  spectrum.  Ideal  bandpass 
filtering  will  recover  our  original  bandpass  waveform.  This 
is  not  true  in  Case  2  of  Figure  16.  Here,  the  upper  cutoff 


Case  2.  v0+<r/2/k<r  and 


,  I 

fa*cr/«  samples/s 


I 


~w. 


Case  3.  v0+<r/2/kcr  and  fa>a/»  samples/s 


frequency  Is  not  an  integer  multiple  of  a  and  If  fs  =  o/n 

samples/s  as  in  Case  1,  the  shifted  spectrums  overlap 
within  the  original  bandpass  spectrum  and  there  is  aliasing. 
Ideal  bandpass  filtering  will  not  recover  our  original 
signal.  A  higher  sampling  rate  is  needed  to  "adjust"  the 
spectra  until  there  is  no  overlapping  spectra  as  in  Case 
3  for  cr/n  <  fs  <  2o/u  samples/s. 

As  stated  before,  the  sampling  frequency  needed  to 
prevent  this  aliasing  will  be  somewhere  between  o/n  and 
2o/u  samples/s.  But  how  do  you  determine  the  actual 
sampling  rate  required?  Refer  to  Figure  17,  assume  that  the 
negative  frequency  spectrum  of  X(w)  is  shifted  just  to  the 
left  of  the  positive  frequency  spectrum  for  some  integer  m, 
the  (m+1)  Integer  shift  must  move  it  completely  to  the  right 
to  avoid  overlap.  In  other  words,  we  require  that 


(-w0  +  o/2)  +  m2nfs 

< 

w0  -  a/2 

(  25  ) 

mws 

< 

2w0-a 

(  26  ) 

ws 

< 

( 2w0 -a ) /m 

(  27  ) 

and 


(-W0  +  o/2)  +  (m+l)2nfs  >  w0  +  a/2 


(  28  ) 


(m+l)ws  >  2w0 


(  29  ) 


ws  >  (2v0)/(m+l) 


(  30) 


Figure  17.  Determining  fs(min)  (6:237-238) 


or 

(2w0)/(m+l)  <  ws  <  (2w0-cr)/m  (31) 

From  this,  ws(min)  can  determined  and  ws  can  be  plotted. 
Let 

M  *  ( (wQ  -  o/2) /o  )  ( 32  ) 

where  (  1  denotes  the  greatest  integer  function.  From 


(  2w0-ct  )  /  [  (  w0-cr/2  )  /a] 


( 


Equation  (34)  is  plotted  in  Figure  18  as  a  function  of  (w 
o/2 )  /o  .  Notice  that  when  (wQ  ♦  a/2)/a  =  integer  , 

sampling  rate  is  always  twice  the  bandwidth,  or  cr/rc 
samples/s.  When  (w0  +  o/2)/o  does  not  equal  an  integer, 
must  be  determined  from  Eq  (34)  or  interpolated  from 
Figure  18. 
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Figure  18.  v8(iain)  •»  a  Function  of  w0  and  a 


Recall  from  our  review  of  low-pass  sampling  that  we 
found  x(t)  in  terms  of  its  samples  in  Eq  (9).  That  equation 
has  the  form 

CD 

x  ( t )  =  E  x(kTs )h(t-kTs )  (9  ) 

k  =  -® 

where  Ts  is  the  period  between  samples,  x(kTs)  are  the 
samples,  and  h ( t )  is  the  low-pass  reconstruction  filter,  or 
the  "interpolating  function." 

It  is  possible  to  extend  Eq  (9)  to  the  case  of  a  direct 
sampled  bandpass  signal.  Since  we  want  to  reconstruct  our 
sampled  signal  by  ideal  bandpass  filtering  over  our  original 
bandpass  spectrum,  we  can  use  the  transfer  pair  in  Figure  19 
and  show  that  our  bandpass  filter  h(t)  is 

h(t)  =  (crTs/Ti)Sa(at/4Ti)cosw0t  (35) 

where  Sa ( x )  =  sinx/x  ,  and  x(t),  in  terms  of  its  uniform 
samples,  becomes 

CO 

x  ( t )  =  (crTs/u)E  x(kTs)Sa[or(t-kTs)/4nJcos[w0(t-kTs)  )  (36) 

k  =  -® 

This  bandpass  sampling  theorem  is  constrained  to  the 
sampling  rate  in  Figure  18  and  Eq  (34). 

Note,  from  Figure  18,  for  narrowband  systems  where 

cr  <<  wQ  =  wQ  +(a/2)  (37) 

it  is  possible  to  sample  at  the  Nyquist  rate  and  still 
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recover  the  bandpass  waveform. 

The  subject  of  direct  bandpass  sampling  followed  very 
closely  the  theory  involved  in  low-pass  sampling.  Our 
bandlimited  bandpass  signal  is  sampled  by  multiplying  it 
with  a  pulse  train.  For  recovery,  a  window  function  (an 
ideal  bandpass  reconstruction  filter)  is  applied  to  our 
properly  sampled  (fs  chosen  to  prevent  aliasing)  signal. 

First-Order  Sampling.  An  alternative  approach  to 
direct  sampling,  also  called  first-order  sampling,  was 
examined  by  John  L.  Brown  (3:613-615).  He  showed  that  the 
first-order  sampling  of  a  bandpass  signal  is  reduced  to  an 
application  of  the  low-pass  sampling  theorem  when  the  upper 
cutoff  frequency  is  an  integer  multiple  of  the  bandwidth. 

He  uses  a  simple  band-embedding  technique  to  restore  the 
positioning  constraint  when  the  upper  cutoff  frequency  is 
not  equal  to  an  integer  multiple  of  the  bandwidth. 

This  proof  uses  only  the  standard  quadrature 
representation  of  the  bandpass  signal  and  the  classical 
sampling  theorem  for  low-pass  signals.  For  real  numbers 
w0,  cr  satisfying  w0  >  a/2  >0  we  define  0(wo,a)  to  consist 
of  all  complex  valued  bandpass  signals  x  which  can  be 
represented  in  the  form 

* 

x(t)  =  (1/2ti)  X  ( w)  exp  (  j  wt )  dw  for  all  t  (38) 

II*  *2 


where 


1 1  =  I1(WQ/Cf)  =  (wo-a/2/wo+cr/2 ) 


(  39) 


1 2  =  -Ii  =  (-w0-ct/2,  -w0+ct/2)  (40) 

and 

j  |  X  ( w )  |  2dw  <  cd  ( 41 ) 

thus  x  possesses  finite  energy  and  all  its  derivatives  are 
continuous.  If  w0+cr/2  =  kcr  for  some  positive  integer  k, 
then  x(t)  can  be  represented  by  samples  taken  at  an  average 
rate  of  ct/ti  samples/s.  If  (w0  +  cr/2)/cr  is  not  an  integer 

(the  general  case),  then  the  bandwidth  a  is  increased  to 
o'  ,  where  ( w0  +  cr/2  ) /cr '  =  (  ( w0+cr/2  ) /ct  ]  ,  and  [  ]  denotes 
the  greatest  integer  function.  For  the  general  case,  an 
average  sampling  rate  between  ct/ti  and  2ct/ti  samples/s  is 
needed  . 

Analysis.  Let  x(t)  be  a  bandpass  signal  having 
all  its  spectral  components  in  the  band  w0-(c/2)  <  |w|  < 
w0+(ct/2)  so  that  (w0,ct)  characterizes  the  position  of  the 
band.  Also,  let  x  f  |3(w0,cr) 

Under  these  assumptions  x(t)  may  be  represented  as 

x(t)  =  p(t)cosw0t  -  g(t)sinw0t  (42) 

for  all  t,  where  p(t)  and  q ( t )  are  low-pass  functions 


bandllmited  to  |w|  <  a/2.  Each  of  these  signals,  p(t)  and 
q(t),  may  be  sampled  at  a  rate  of  a/2n  samples/s. 

This  analysis  is  broken  into  two  phases.  In  the  first 
phase  we  assume  w0  +  a/2  =  ka  ,  where  k  is  a  fixed  positive 
integer  and  in  the  second  phase  (the  general  case)  we  assume 
wQ  >  a/2  >  0  (arbitrary  band  positioning). 

Case  i.  Where  x(t)  is  defined  by  Eq  (42)  and 

x(nT)  =  p ( nT )cos ( nw0T )  -  q ( nT ) s i n ( nw0T )  (43) 

Assume  T  =  n/a  corresponding  to  a  sampling  rate  of  cr/n 
samples/s.  So, 

x(nT)  =  p  (  nT  )  cos  (  nrtw0/a  )  -  q  (  nT  )  s  i  n  (  nnw0 /a  )  (44) 

and  with  w0  +  a/2  =  ka  ,  where  k  is  fixed  positive 
integer  (w0/a  =  (k-1/2)  =  ( ( 2k -1 ) /2 ) } 

x(nT)  =  p ( nT ) cos ( nn ( 2k -1 ) /2 3 

-  q ( nT ) s in [ nn ( 2k -1 ) /2 )  (45) 

Where  n  is  even,  n  =  2v  ,  and  letting  =  2T  =  2 n/a 

the  sine  term  goes  to  zero  and 

x ( 2 vT )  =  x(2vT1/2)  =  x ( vTx )  =  p(vT!)(-l)V  (46) 

Where  n  is  odd,  n  =  2v-l  ,  the  cosine  term  goes  to  zero 
and 


x( ( 2v-l ) T 1  =  x ( (2v-l)T1/2) 


(47) 


x(t)  =  E  p(2vT)$  (t-2vT)cosw  t 


v= -® 


E  q(  2vT-T)*0-(t-2vT+T)sinv0t 

V  =  -oo 


but  from  Eq  ( 46 ) 


p ( 2vT )  =  x { 2vT ) ( -1 ) 


and  from  Eq  ( 48 ) 


q ( 2vT-T )  =  xl(2v-l)T](-l) 


v+k 


Plugging  these  into  Eq  (53) 


x ( t )  =  E  x(2vT)(-l)  ♦a(t-2vT)cosw0t 

V  =  -00 


v+k. 


-  E  x[  (2v-l)tJ  (-1)  ’'♦CT(t-2vT+T)sinv0t 


v = -  ® 


Now,  letting  v  =  n 


x(t)  =  E  x(2nT) (-l)nta(t-2nT)cosw0t 
n=-« 


-  E  x( (2n-l)t J (-l)n+k*CT(t-2nT+T)sinw0t 
n=-« 


Since  w0t 


=  ( 2k-l )n/2 


cosw0t  =  cosw0(t-2nT) 


and  since 


( -1  )n  +  *Csinw0t  =  cosw0  [  -  (  2n- 1 } T  ] 

Eq  (57)  becomes 

00 

x(t)  =  E  x(vT)*a(t-vT)cosw0(t-vT) 

-OD 

which  is  the  desired  form  of  the  expression  for  the  case 
wc  +  cr/2  =  kcr 

Case  i i .  For  the  general  case,  we  assume 
arbitrary  band  position  and  w0  >  cr/2  >  0.  Define 

r  =  [ ( w0+a/2 ) /a ] 

with 

r  <  ( w0+ct/2 ) /a  <  r+1 

Holding  the  upper  cutoff  frequency,  wD  +  cr/2  , 


constant  (Figure  20)  and  increasing  the  bandwidth  a  to 

a'  such  that  (w0  +  a/2)/a'  =  r 

Now,  !‘j 

w0'  =  (w0+cr/2)-cr'/2 

(  6  3  )  V 

1 1  ( «o  '  /  °  '  )  >  1 1  ( w0  ,  CT  ) 

(64) 

Ci 

1 2 ( wo ' • °  '  }  >  1 2 ( wo ' CT  } 

(65) 

and  x  e  0 (wQ ' ,  cr ' ) 


!ll(w0/a) 


Il(v0t,o'  ) 


Hold  the  upper  cutoff 
frequency  fixed  while 
increasing  the  bandwidth 
from  o  to  o' 


Figure  20.  Embedding  Technique  for  First-Order  Sampling 

The  upper  cutoff  frequency  w0  ♦  o/2  is  a  multiple  of 
o’  and  o,  T,  and  wQ  are  replaced  in  Eq  (60)  by  o'  , 

T'  *  x/o '  ,  and  wQ '  respectively. 

We  now  have 

x(t)  =  x(nT’ )♦_, (t-nT* )cosw  ' (t-nT’ )  (66 


So  a  representation  of  x(t)  Is  possible  with  uniform  samples 
at  a  rate  of 

fs  =  1/T '  =  cr'/m  =  ttor/[aln  (67) 

samples/s  where  a  =  ( w0  +  a/2)/a 

Summary.  Uniform  sampling  of  a  bandpass  signal  at  an 
average  rate  of  a/n  samples/s  will  suffice  to  determine 
the  signal  for  all  values  of  time  only  when  the  upper  cutoff 
frequency  is  an  Integer  multiple  of  the  bandwidth.  For  the 
case  where  the  upper  cutoff  frequency  is  not  an  integer 
multiple  of  the  bandwidth,  an  average  rate  between  a/n  and 
2a/n  samples/s  is  required.  For  direct  sampling,  this  rate 
can  be  determined  from  Eq  (34)  or  from  Figure  18.  For 
first-order  sampling.  Brown  demonstrated  a  simple  band 
embedding  technique  which  can  be  used  to  determine  the 
minimum  sampling  rate  for  the  signal.  Brown's  approach  used 
only  the  quadrature  representation  of  the  bandpass  signal 
and  the  classical  Shannon  sampling  theorem  in  his  proofs. 

The  sampling  rate,  fs,  for  the  general  case  is  given  by 
Eq  (67). 

pth-order  Sampling 

I nt r oduct  1  on .  Kohlenberg  extended  the  idea  of  first- 
order  sampling  to  pth-order  sampling  (8:1432-1436).  Unlike 
first-order  sampling  which  involves  a  single  train  of 
uniformly  separated  samples  of  the  bandpass  waveform,  pth- 
order  sampling  Involves  p  trains  of  samples.  The  trains  are 


slightly  displaced  In  time  but  occur  at  the  same  frequency 
(10:40) . 


The  following  development  discusses  pth-order  sampling 
followed  by  second-order  sampling  (the  case  where  p  =  2). 

To  date,  only  quadrature  sampling,  a  special  case  of  second- 
order  sampling,  has  been  shown  to  be  a  practical  approach  to 
pth-order  bandpass  sampling.  For  that  reason,  quadrature 
sampling  will  be  covered  in  greater  detail. 

The  second-order  development  will  lead  to  a 
representation  of  x(t)  in  terms  of  its  uniform  samples  just 
as  in  the  previous  developments;  however,  in  this  case,  x(t) 
will  be  represented  by  two  trains  of  uniform  samples. 

pth-Order  Sampling  Analysis.  (8:1442-1443)  Kohlenberg 
considered  the  spectra  of  amplitude  modulated  sequences  of 
pulses  of  the  form 


g ( t )  =  Z  x(an)s( t-an ) 
n  =  -® 


(68) 


where 

s(t)  gives  the  shape  of  the  basic  pulse,  and 
x(t)  represents  the  modulating  amplitude 

Here,  the  function  g(t)  represents  the  first-order  sampling. 
Kohlenberg  defined  pth-order  sampling  as  a  function 


P 

g ( t )  =  z  g  A ( t )  (69) 

i  =  l 
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( 


p 

g(t)  =  E  IE  x(a1n+k1 )si (t-a1n-k1  )  1 
i = 1  n=-« 

where  the  ith  sampling  of  g^ft)  has  step  size  a^,  phase 
kj,  and  sampling  function  S}(t). 

To  find  the  Fourier  transform  of  g(t)  consider 


g^t)  =  E  x(a1n+ki  )si  (t-ain-kl ) 
n=-® 


Using  the  Dirac,  6,  function  this  becomes 


=  j 


+  00 


gA(t)  =  E  xtYjsj (t-Y)d(Y-ain-ki )dY 


Interchanging  the  order  of  summation  and  integration 


•  r: 


x ( Y ) s j ( t-Y ) {E  d ( Y-a i n-k  j  ) ) dY 

n  =  -« 


This  summation  in  brackets  represents  a  periodic  function 
with  period  a*.  Using  its  Four ier-ser ies  expansion 


-1 


E  a(Y-ajn-ki)  =  E  a^  exp{ j2nn(Y-ki )/aj) 
n  =  -oo  n=-® 


and  rewriting 


9l<t) 


I 


+  00 


-1, 


x(Y)sj(t-Y)ai  {E  exp{  j2tin  ( Y-k  j  )/a  j  )dY  ( 
.  n  =  -« 
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gi(t)=ai  E  {  s  i  ( t-Y  )  x  ( Y  )  exp  (  j  2unY/a  ^  )  dYexp  (  -  j  2Tink  i /a  ^  )  } 

n  =  -®  J-o, 

(76  ) 

Taking  the  Fourier  transform  (F.T.)  of  both  sides  of  this 
equation  and  using  the  fact  that 

F.  T . { x ( Y ) exp ( j 2nnY/a  i  )  }  =  X(w+2nn/ai)  (77) 

-l  00 

G^tw)  =  a^  ~E  exp ( - j 2mnk ^ /a i ) 3 j ( w) X ( w+2nn/ai )  (78) 

n  =  -co 

Summing  over  all  "i"  yields 

P  00 

G ( w )  =  E  a^  S  ^ ( w ) { E  X ( w+ 2mn/a ^ ) exp ( - j 2unk j /a ^ ) }  (79) 

1  =  1  n  =  -  <*> 

This  pth-order  sampling  theorem  gives  the  Fourier  transform 
of  g(t)  in  terms  of  the  sampled  function  x(t)  and  the 
sampling  parameters. 

Second-Order  Sampling  Analysis.  (13:907)  Second-order 
sampling  is  a  specialized  case  of  pth-order  sampling  (p=2). 
Two  equally  effective  implementations  of  second-order 
sampling  are  pictured  in  Figure  21.  In  one  case  the 
bandpass  signal  is  passed  through  two  channels,  one  channel 
is  delayed  in  time  and  both  channels  are  then  simultaneously 
sampled.  In  the  other  implementation,  the  bandpass  signal 
is  once  again  fed  through  two  channels;  however,  this  time 
one  of  the  sample  trains  is  delayed  (Figure  22).  The 
samples  occur  at  the  same  frequency  but  are  delayed  in  time 


from  each  other  . 


Sample  i  rain  2 


•  2  T,  ♦  r  Or  IT  ♦  r  AT.  ♦  r  r>Tt  +  r  8  T,  ♦  r  / 
ibi 


Figure  22.  Sample  Trains  For  Second-Order  Sampling  (10:41) 

As  mentioned,  second-order  sampling  yields  two  sets  of 
uniformly  spaced  samples  which  are  interleaved,  i.e., 

x  ^ ( mT )  =  x ( mT )  (80) 

and 

x  2  (  mT  )  =  x(mT-ct)  (81) 

where  T  is  the  temporal  offset  of  the  staggered  sequence. 

The  sequence  x^fmT)  has  the  Fourier  transform  given  by 


X]_(v)  =  X  ( w ) 


(82) 


and  the  transform  of  x2(mT)  is 

00 

X2(w)  =  (1/T)E  X( w+2nm/T)exp{ - j [w+ (2nm/T) ]a}  ( 

m=  -  ® 

The  use  of  x2(mT)  introduces  a  new  degree  of  freedom  over 
first-order  sampling  which  can  be  used  to  eliminate 
al iasing . 

The  original  waveform  can  be  recovered  by  properly 
filtering  and  combining  the  two  sequences,  i.e., 

OD  00 

x(t)  =  £  x  jl  (  mT  )  hp  ( t-mT  )  4-  c  x2  (  mT  )  hp  ( -t  +  mT-a  )  ( 

m=~®  m  =  -® 


For  exact  reconstruction  an  ideal  bandpass  filter  is 
required  with  frequency  transfer  characteristics,  HD(w) 
given  by 


Hd(w) 


where 


0,  0  <  w  <  w0  -  ( n/t ) 

Ai,  w0-(n/T)  <  w  <  {{2n[  ( 2k - 1 ) / 2 ] ) /T } -w0 
A2,  { 2 [ ( 2k-l )/2 1 /T } -w0  <  w  <  w0+(n/T) 

0,  w0  +  ( jt/T )  <  w  <  2m/T  ( 


A:  =  { exp  {  -  j  [  (  (  2k -1 ) /2  )  (3 -ts/2  3  }/(2/T)sin{  [  (  2  k  - 1 )  /  2  1 13  J  ( 


A2  =  { exp  [  -  j  (  kfl-Ti/2  )  ]  }  /  (  2/T )  s  1  n  (  k!3  ) 


(3  =  2na/T 


( 


and  k  is  the  largest  integer  for  which 


r, 

r. 

f: 


C 


[(k-l)/T]  <  f0-(l/2T)  (89) 

Conventional  Quadrature  Sampling 

I ntroduct  i  on .  As  mentioned  in  Chapter  I,  it  is 
possible  to  preprocess  a  signal  prior  to  sampling  to  insure 
being  able  to  sample  at  the  Nyquist  rate.  Recall,  that  a 
bandpass  signal  can  be  represented  in  the  form 

x(t)  =  p(t)cosw0t  -  q(t)sinw0t  (42) 

where  p ( t )  and  q(t),  the  in-phase  and  quadrature  components, 
respectively,  are  low-pass  signals  with  spectra  occupying 
the  symmetric  band  jwj  <  a/2  ;  effectively,  the  bandpass 

signal  is  specified  by  a  pair  of  low-pass  signals. 

The  objective  of  conventional  quadrature  sampling  is  to 
recover  the  low-pass  components,  p ( t )  and  q(t),  and  to 
sample  them  at  the  Nyquist  rate.  Because  of  the 
preprocessing  of  x(t),  the  low-pass  components  are 
independent  and  may  be  sampled  independently  (10:36).  This 
means  that  there  is  only  the  Nyquist  limit  to  satisfy.  We 
no  longer  have  to  be  concerned  with  two  bandpass  spectrums 
overlapping  or  the  relationship  between  the  upper  cutoff 
frequency  and  the  bandwidth.  Samples  of  the  low-pass 
components  will  uniquely  determine  our  original 
bandpass  signal. 

Analysis.  A  common  approach  to  obtaining  the 
quadrature  components  is  to  employ  quadrature  demodulation 


-  a  -~n  Jrn  -'p  m\  m*l 


44 


as  shown  in  Figure  23  (13:907-908).  In  one  channel  x(t)  is 
modulated  by  cosw0t,  or 


x(t)cosw0t  =  [p(t)cosw0t  -  q(t)sinw0t)cosw0t  (90 

2 

x(t)cosw0t  =  p(t)cos  w0t  -  q ( t ) s inw0tcosw0t  (91 

Recall, 

2 

cos  A  =  1/2  +  (l/2)cos2A  (92 

sinAcosB  =  ( 1/2 ) ( s in ( A+B )  +  sin(A-B)]  (93 

Here  A  =  B,  so 

sinAcosA  =  (l/2)sin2A  (94 

x(t)coswQt  =  (1/2) tp(t)  +  p(t)cos2w0t  -  q(t)sin2w0t)  (95 

The  filter  removes  the  2w0  terms,  and 

x(t )cosw0t  «  p(t)/2  (96 

Similarly,  q(t)  is  obtained  from  the  sine  channel. 

2 

x(t)sinw0t  =  p ( t ) s i nw0tcosw0t  -  q(t)sin  wQt  (97 

sin2w0t  =  ( 1/2 ) [ l-cos2w0t ]  (98 

x(t)sinw0t  =  ( 1/2 ) [ p ( t ) s in2w0t  -  q(t)  +  q(t)cos2wct)  (99 
Once  again,  filtering  out  the  2w0  terms 

x(t)sinw0t  ~  -q(t)/2  (10C 
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Both  channels  are  then  sampled  at  a  rate 
of  cr/2n  samples/s.  The  samples  are  written  as 

CD 

p(t)  =  E  p(kTs)<5(t-kTs)  (101) 

k  =  -co 

CD 

q  ( t )  =  E  q  ( kTs  )  <5  ( t-kTs  )  (102) 

k  =  -® 

These  equations  have  the  familiar  form  of  a  low-pass  signal 
in  terms  of  its  samples.  The  next  step  in  the  process  is  to 
apply  the  samples  to  the  reconstr uct i on  filters  as  in  Figure 
23b.  These  reconstr uct i on  filters  are  the  spectral  windows 
referenced  earlier.  Signals  p(t)  and  q(t)  are  modulated 
once  again  to  yield  the  familiar  form  of  Eq  (42) 

x(t)  =  p(t)coswQt  -  q(t)sinw0t  (42) 

Summary.  Conventional  quadrature  sampling  allows  us  to 
sample  an  arbitrary  bandlimited  bandpass  signal  at  a  total 
rate  of  ct/ti  samples/s,  using  preprocessing,  regardless  of 
the  ratio  of  the  center  frequency  and  the  bandwidth,  and 
recover  x ( t ) . 

Quadrature  Sampling 

Intr oduct 1  on .  This  approach  to  quadrature  sampling  is 
illustrated  in  Figure  24.  Instead  of  preprocessing  to 
determine  the  low-pass  components,  p ( t )  and  q(t)  are 
determined  from  samples  of  the  bandpass  signal  and  its 
quarter-wavelength  extension  (7:1453-1454). 


Figure  24.  Quadrature  Sampling 

Three  cases  will  be  discussed.  In  Case  i  it  is  assumed 
that  2w0/cr  =  k  ,  where  k  is  an  integer  t  1.  It  will  be 
shown  that  sampling  at  an  average  rate  of  a/u  samples  is 
sufficient  to  recover  the  bandpass  signal.  Case  ii  assumes 
only  that  0  <  a/2  i  wQ  <  «  ;  i.e.,  2w0/a  may  or  may 
not  equal  k.  A  technique  is  applied  where  Ii(wQ,a)  is 
embedded  in  a  larger  band,  Ii(v0,a'),  having  the  same 
center  frequency  w0,  and  a  bandwidth  a'  i.  a  .  In  Case 
iii,  I i ( wQ , a )  is  embedded  in  an  "optimal"  band,  Ii(w0',o"), 
that  requires  2w0’/a"  to  be  an  integer.  Cases  ii  and  iii 
allow  sampling  at  an  average  rate  between  a/n  and  2 a/n 

samples/s  as  will  be  explained  in  detail  in  the  following 
paragraphs . 

Case  1.  (2:1660)  Quadrature  sampling's  goal  is 
to  recover  the  low-pass  components  p(t)  and  q(t)  directly 
from  samples  taken  of  the  bandpass  signal  x(t)  and  its 


quarter-wavelength  translation  x(t-*/2w0)  (2:1659).  Case  1 
assumes  the  constraint  2w 0/cr  =  k  ,  where  k  Is  an 
integer  t  1. 

Recalling  Eq  (42) 

x(t)  =  p(t)cosw0t  -  q ( t ) s inw0t  (42) 

In  terms  of  its  samples  x(t)  becomes 

xfnTj)  =  pfnTjJcosfnwQTi)  -  q  ( n^ )  s  in  ( nw0T! )  (103) 

Letting  =  2re/cr  and  2w0/cr  =  k 

x(nTi)  =  p(nTi)cos(n(ak/2)  (2n/cr)> 

-  q(nT1)sin(n(ok/2) (2tc/ct) }  (104) 

x(nT1)  =  p ( nT^  ) cos  ( nxk  )  -  q ( n^  ) sin  ( nnk  )  (105) 

or, 

x(nTx)  =  ( -1 )knp(nT1 )  (106) 

x(nT1-n/2w0)  =  (-l)knq(nT1-n/2w0)  (107) 

These  equations  show  that  samples  of  p(t)  at  a  rate  of 
1/Tl  =  cr/2u  samples/s  are  available  from  the  original 
bandpass  waveform  x(t),  as  are  samples  of  q(t-x/2w0)  at  o/2n 
samples/s.  Recall  that  p(t)  and  q(t)  are  bandlimited  to 
| w J  <  o/2  ,  which  implies  that  they  can  be  reconstructed 

from  their  samples  at  the  low-pass  rate. 

From  the  Shannon  low-pass  sampling  theorem 


(108) 


p(t)  =  E  p(nTi )*(t-nTi ) 
n=-® 

where 

♦(t)  =  (sin(at/2) )/(<rt/2) 

So,  plugging  Eq  (106)  into  Eq  (108)  yields 

<*>  , 

p(t)  =  E  (-l)Knx(nT1)*(t-nT1) 
n  =  -® 

Similarly, 

00 

q(t-a)  =  E  q(nT1-a)*(t-nT^) 
n=-® 


where  a  is  real. 

Let  t  =  t  +  a 

00 

q(t)  =  E  q(nT^-a)*(t+a-nTi) 
n= -® 

For  a  =  tc/2wq 


q ( t )  =  E  q(nT1-n/2w0)e{t+(n/2w0)-nT1} 
n=-« 

Plugging  Eq  (107)  into  Eq  (112) 

ao 

q(t)  =  E  (-l)knx(nT1-n/2w0)*{t+(n/2w0)-nT1) 
n=-® 


(109  ) 


(110) 


(111) 


(112) 


(113) 


(114) 
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So,  by  plugging  Eqs  (110)  and  (113)  Into  Eq  (42) 


x(t)  =  E  ( -1  )  knx  (  nTi )  <t>  ( t -nT]_ )  coswQt 

n  =  -® 

00  kn 

E  ( -1 )  x  (  nTi -ti/2w0  )  ♦  { t  +  ( u/2w0  ) -nTi }  s  inwQt  (115) 

n  =  -® 

which  is  the  desired  quadrature  sampling  formula  for  the 
case  2w0/a  =  k  =  integer  >  1 

The  low-pass  components  p(t)  and  q(t)  were  sampled  at  a 
rate  of  cr/2u  samples/s  yielding  an  average  rate  of  ct/ti 
samples/s  for  the  bandpass  signal. 

Case  i  i  .  (2:1660)  For  this  case  assume  only  that 

®  >  w0  >  cr/2  >  0  .  In  other  words,  2w0/cr  need  not  always 

be  an  integer.  Here  Ii(w0,cr)  is  embedded  in  a  larger  band 
II ( w0,a ' )  which  has  the  same  center  frequency  and  a 
bandwidth  cr '  >  a  as  shown  in  Figure  25.  Holding  w0 
constant,  and  increasing  a  (if  necessary)  to  the  value  o' 
such  that  2w0/a'  will  be  an  integer,  yields 

2w0/cr'  =  [2w0/o]  =  s  =  integer  (116) 


where 


s  <  2w 0/o  <  s+1  (117) 

The  results  of  Case  i  now  apply  with  a  ,  T]_  ,  k  , 

replaced  with  a'  ,  '  =  2n/cr'  ,  s  respectively.  So, 


Il<  Vo/0r) 


I— cr-H 


> - 1 - «*w 


Hold  the  center 
frequency  fixed  while 
increasing  the  bandwidth 
from  o  to  o' 


•  r« 


( -1 )  x ( nT: • )*1(t-nT1 ' ) 


(1 


p(t)  =  Z 

n  =  -  os 


q(t)  =  Z  (-l)Snx(nT1'-Ti/2w0)4>1{t+(7i/2w0)-nT1' }  (1 

n  =  -® 

where 

♦  l_(t)  5  {sin(a't/2)  }/(cr't/2)  (1 

Sampling  x(t)  and  x(t-n/2w0)  at  the  rate 

1/T1'  =  o'/ 2k  =  {  (2w0/cr)/I2w0/cr]  }  (a/2n)  (1 

or 

1/TX '  =  (2a/[2a]J (a/2u)  (1 

yields  an  average  sampling  rate 

(2«/t  2a]  )  (cr/Ti)  (1 

samples/s  for  the  reconstruction  of  x(t)  from  Eq  (42) 

x(t)  =  p(t)cosw0t  -  q(t)sinw0t  ( 

When  2w Q/a  is  an  integer  Eq  (123)  becomes  cr/u  samples/s, 
which  is  the  equivalent  to  Case  1  as  expected. 

Case  i i 1 .  (2:1661)  In  this  case  it  is  not 
necessary  to  hold  w0  or  a  fixed.  Here,  x(t)  is  a  general 
bandpass  signal,  w0  doesn't  have  a  strict  sense  of  center 
(angular)  frequency,  but  only  a  sense  of  "average  of 


extremes" 


(w0-ct/2  and  w0+a/2). 


Again,  assume  that  •  >  w0  >  a/2  >  0  ,  and  that  2w 0/a 

may  or  may  not  equal  an  integer.  The  objective  is  to 
determine  a  generalized  result  by  allowing  another  degree  of 
freedom  in  the  embedding  technique  of  Case  11. 

For  this  embedding  technique,  the  upper  cutoff 
frequency,  w0+a/2,  is  held  invariant  while  increasing  the 
bandwidth  (Figure  26).  An  "optimal"  covering  band 
ll(wo,,cr">  ls  determined  where  it  is  required  that  2w0'/cr"  = 
integer.  The  analysis  follows. 

Let  s  =  [2w0/cr]  be  the  greatest  integer  contained  in 
2w0/cr  and  s  <  2w0/cr  <  s  +  1  .  From  Figure  26,  it  is 

obvious  that 

wQ’  *  (w0+o/2)-or"/2  (124) 

where  a"  >  a 

To  determine  a "  ,  set 

2w0'/cr"  *  2{w0+(cr/2)-<7"/2}/a"  =  s  =  12w0/ct)  (125) 


or 


2w0'/a"  =  ( 2w0+a-cr"  ) /a"  =  s  =  [2w0/ol  (126) 

Solving  for  a" 

a"s  =  2w0+o-cr"  (127) 

o"s+a"  =  2w0+cr  (128) 


a" (s+1)  =  2w n+a 


(129) 


1 1  (  v0 ,  cr ) 


Hold  the  upper  cutoff 
frequency  fixed  while 
increasing  the  bandwidth 
from  a  to  a" 


Figure  26.  Embedding  Technique  -  Case  iii 


■»  s  *. 


o"  =  {cr(2w0/<r)  +  l}/U2v0/cr]+l}  (13 

Now,  let  Y  =  w0/cr 

a"  =  (a ( 2Y  +  1 ) } /[ 2Y  +  1 )  (13 

If  Y  >  1  then  cr"  <  2a  and  w0-a/2  <  w0'  <  w0. 
Further,  v0'+o"/2  =  w0+o/2  ;  l.e.,  the  upper 
cutoff  frequency  of  the  covering  band  remains  Invariant. 

The  low-pass  components  p(t)  and  q(t)  with  respect  to 
the  center  frequency  w0'  have  the  form 


p(t)  = 


00 

£  (-1) 
n  =  -® 


sn 


xinTi" )*2 (t-nTi” ) 


(13 


q  ( t )  = 


00 

E  (-1) 
n  =  -« 


an 


x(nT1"-n/2w0*  )  *2  i  t+  U/2w0  '  J-nTj,"}  (13 


where 


♦2(t)  =  {sln(a"t/2) }/(a"t/2) 


(13 


and 


Tj,"  =  2Ti/a" 


(13 


With 


x(t)  =  p(t)cosw0’t  -  q(t)sinw0’t 


(13 


»  • 
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>  \  .V 


x  ( t )  =  Z  ( -1 ) snx ( nTi" } $2 (t -nTj " ) cosw0 ' t 
n  =  -«o 

as 

Z  ( -1  )  snx  (  nTi  " -ti/2w0  ’  )  $2  ft  +  (  n/2w0  '  ) -nT]_  "  }  s  inw0  '  t 

n  =  -od 

(138) 

the  average  sampling  rate  in  this  reconstruction  is 

{  (2Y  +  D/I2Y  +  1D  (o/n)  (139) 

samples/s .  Si  nee 

( 2Y  +  1 ) / [ 2Y  +  1  ]  <  2Y/(2Y]  (140) 

the  sampling  rate  is  improved  over  that  in  Case  ii. 

Summary.  In  Case  i  an  average  sampling  rate  of  o/n 
samples/s  was  determined  for  the  reconstruction  of  x(t)  from 

x(t)  =  p(t)coswQt  -  q(t)sinwDt  (42) 

for  the  situation  where  the  upper  cutoff  frequency  was  an 
integer  multiple  of  the  bandwidth.  When  this  constraint 
doesn't  hold,  the  average  sampling  rate  was 

{ ( 2Y  ) / ( 2Y ] } ( o/n  )  samples/s  as  determined  in  Case  ii.  Since 
1  <  2y/[2y]  <  2  for  all  Y  >  1/2  ,  the  average  rate 

necessary  for  reconstruction  of  x(t)  varies  from  o/n  to  o/2n 
samples/s.  Lastly,  when  the  frequency  wQ  is  allowed  to 

vary,  a  further  improvement  in  the  sampling  rate  can  be 
obtained,  the  average  rate  being  given  by  Eq  (139). 


Complex  Sampling 


Intr oduct 1  on .  Another  possibility  for  maintaining  the 
(average)  sampling  rate  equal  to  the  minimum  value  of  a/n 
samples/s  is  to  sample  both  x(t)  and  its  Hilbert  transform 
at  a  rate  of  1/TS  =  a/2m  samples/s  (10:37).  Since  x(t)  and 
its  quadrature  function  are  Hilbert  transforms  of  each 
other,  the  low-pass  components,  p(t)  and  q(t),  can  once 
again  be  determined,  this  time  from  samples  of  x(t)  and  its 
Hilbert  transform,  x^t). 

Ana  I ys 1 s .  (13:906  )  This  technique,  illustrated  in 

Figure  27,  is  equivalent  to  forming  the  analytic  signal 


Zx(t)  =  x ( t )  +  jxH(t) 


1 1 41 ) 


x  ( t )  - f 


•  x ( nT_ ) 


p(t)  =  £  6( t-nTs) 


xH(t) 


Hilbert 

Transform 


XH ( nTs ) 


Figure  27.  Complex  Sampling  (13:905) 


The  Hilbert  transform  of  x(t),  xu(t),  is  obtained  by  the 


following  operation 


xH(  t )  =  - jsgn ( w) X( w) 


(142) 


The  Fourier  transform  of  the  analytic  signal,  Zx(t),  is 


2v(t)  = 


2X ( w ) , 
X  ( w  )  , 
0, 


w  >  0 
w  =  0 
w  <  0 


(143  ) 


For  the  bandpass  case,  X(0)  =  0  and  Zx(t)  =  2X(w)U(w),  where 
U  is  the  step  function. 

In  terms  of  the  positive  and  negative  frequency 
spectra,  X+(w)  and  X_(w),  respectively,  which  are 
illustrated  in  Figure  28, 


X ( w )  =  X+(w)  +  X_(w) 

(144  ) 

where 

X  (  w )  ,  w  >  0 

X  +  (w)  = 

_  0,  otherwise 

X  (  w  )  ,  w  <  0 

X_(v)  = 

0,  otherwise 

(145) 

(146  ) 

Thus , 

Zx(w)  =  2X+(w) 

(147) 

The  analytic  signal 

Zx(t)  =  x(t)  +  jxH ( t  ) 


(148  ) 


has  a  spectrum  confined  to  the  positive  frequencies 
wQ  -  a/2  <  w  <  wQ  +  a/2.  So  we  can  write 

CD 

Zx(t)  =  E  Zx ( nTs ) s i nc [ a ( t-nTs ) /2 3  exp [ j w0 ( t -nTs )  ]  (149) 

n=  -  ® 

and  from  this, 

x ( t )  =  Re { Zx ( t )  }  (150) 

OD 

x(t)  =  E  { x ( nTs ) s inc (a ( t-nTs ) /2 3  cos [ wQ ( t-nTs ) 3 

n  =  -oo 

-  x^ ( nTs ) s inc [ a ( t-nTs ) /2 3 s in [ w0 ( t-nTs ) 3 }  (151) 

Once  again  we  have  the  familiar  form  of  Eq  (42) 

x(t)  =  p(t)cosw0t  -  q(t)sinw0t  (42) 

Summary.  Obviously  there  is  not  a  strict  separate 
interpolation  in  Eq  (151)  since  the  "carriers"  are 
cos [ w0 ( t-nTs ) 3  and  s in [ wQ ( t-nTs ) 3 .  But,  once  again,  the 
bandpass  signal  is  represented  in  terms  of  its  uniform 
samples . 

Chapter  Summary 

As  mentioned  in  the  Introduction,  the  purpose  of  this 
chapter  was  to  develop,  theoretically,  the  bandpass  sampling 
techniques.  In  each  case,  x(t)  was  developed  in  terms  of 
its  uniform  samples.  In  addition,  the  minimum  sampling  rate 
was  developed  for  each  technique. 


Ill . 


Results 


I ntroduct Ion 

One  "figure  of  merit"  for  comparing  bandpass  sampling 
techniques  is  to  compare  the  minimum  sampling  rate  possible 
for  each  case.  The  following  is  a  review  of  the  key  results 
o  f  Chapter  1 1  . 

Direct  sampling  can  be  accomplished  at  a  minimum  rate 
of  fs  =  cr/rc  samples/s  only  when  w0  +  a/2  =  kcr  (k  = 
integer).  Otherwise,  the  minimum  sampling  rate  is  governed 
by  the  exact  relationship  between  the  bandwidth  and  the 
center  frequency  and  can  be  determined  from  Eq  (34)  or  from 
Figure  18  . 

For  conventional  quadrature  sampling,  the  low-pass 
components  can  each  be  sampled  at  a  minimum  rate  of 
fs  =  a /2ti  samples/s  for  a  total  (average)  rate  of  fs  =  a/n 
samples/s.  This  approach  allows  arbitrary  band  positioning. 
In  other  words,  the  relationship  between  the  center 
frequency  and  the  bandwidth  is  not  considered  when 
determining  the  sampling  rate.  The  preprocessing  allows  the 
low-pass  components  to  be  independently  sampled  according  to 
the  low-pass  sampling  theorem  (10:36).  Since  we  are  usually 
interested  in  the  amplitude,  A(t),  and  the  phase,  #(t),  of 
x(t),  the  time-varying  nature  of  our  signal  requires  that 
samples  of  p(t)  and  q(t)  be  taken  at  the  same  instant 


(18:731)  . 

Quadrature  sampling  allows  sampling  at  the  average  rate 


of  fs  =  cr/n  samples/s  when  2wQ  =  kcr  (k  =  Integer). 

Otherwise  fs  =  { (  2Y  +  1 ) / 1  2Y  +  1  ] )  (cr/m )  samples/s,  where  Y  = 
w 0/a  .  As  in  direct  sampling,  the  minimum  sampling  rate  is 

a  function  of  the  relationship  between  the  center  frequency 
and  the  bandwidth. 

Another  possibility  for  maintaining  the  (average) 
sampling  rate  equal  to  its  minimum  value  or/tr  samples/s  in 
any  case  (arbitrary  band  position)  is  to  sample  both  x(t) 
and  its  Hilbert  transform,  xH(t),  at  a  rate  of  fs  =  a/2it 
samples/s . 

The  last  three  approaches  offer  a  sampling  rate 
reduction  of  up  to  one-half  over  direct  sampling,  but  none 
of  the  approaches  offer  a  significant  rate  reduction  over 
each  other.  So,  how  do  you  chose  a  sampling  scheme?  The 
best  answer  to  that  question,  ultimately,  is  determined  by 
how  the  sampling  scheme  will  fit  into  the  overall  system 
(receiver)  design.  There  are  many  factors  that  must  be 
considered  when  choosing  a  sampling  scheme.  For  instance, 
if  your  design  is  being  driven  by  cost,  it  may  be  cheaper  to 
use  a  single  A/D  and  use  direct  sampling.  Or,  using  "off 
the  shelf"  equipment  may  keep  the  cost  down  and  at  the  same 
time  drive  your  design  to  match  the  available  hardware.  If 
the  system  design  requires  down-conversion,  direct  sampling 
or  the  conventional  quadrature  approach  may  be  appropriate. 

Since  there  is  no  real  "figure  of  merit"  to  compare  the 
individual  sampling  techniques  (unless  you  base  your 


decision  on  the  minimum  sampling  rate  possible),  selection 
of  a  sampling  scheme  should  be  made  only  after  an 
examination  of  what  it  actually  takes  to  implement  these 
techniques  in  hardware. 

A  quick  glance  at  the  figures  in  Chapter  I  gives  a 
pretty  good  idea  of  the  hardware  components  needed,  but 
there  is  more  to  consider.  Recall  in  low-pass  sampling,  the 
theoretical  minimum  sampling  rate  was  shown  to  be  twice  the 
highest  frequency.  That  may  be  true  in  theory,  but  in 
practical  applications  the  rate  is  higher.  A  good  rule  of 
thumb  has  been  shown  to  be  five  times  the  highest  frequency 
(14:53).  The  reason  for  the  difference  is  that  in  the 
theoretical  derivation  all  things  were  considered  ideal. 
Signals  had  bandlimited  spectrums  and  filters  could  pass 
or  reject  whatever  frequencies  were  required.  But  the  fact 
is,  signals  are  not  ideally  bandlimited,  and  since  bandpass 
sampling  rates  are  typically  determined  based  upon 
consideration  of  the  bandwidth,  some  thought  needs  to  be 
given  to  the  definition  of  bandwidth.  A  few  of  the  more 
popular  definitions  are  given  below  (4:61-62). 

Half-power  bandwidth  is  the  interval  between  the 
frequencies  where  the  signal  has  dropped  to  half  power  (or 
3  dB  below  the  peak  value). 

The  equivalent  rectangular  bandwidth  is  defined  as  the 
bandwidth  which  satisfies  the  relationship  P  =  WnS(fc), 
where  P  is  the  total  signal  power  over  all  frequencies,  Wn 
is  the  equivalent  noise  bandwidth,  and  S(£c)  is  the  value  of 


the  signal  at  the  band  center  frequency  (assumed  to  be  the 
maximum  value  over  all  frequencies). 


The  null-to-null  bandwidth  is  the  width  of  the  main 
spectral  lobe,  where  most  of  the  power  is  contained. 

There  are  other  definitions  of  bandwidth  available,  but 
the  point  is,  there  are  many  definitions  of  bandwidth  and 
the  definition  you  choose  has  a  direct  affect  on  the 
sampling  rate,  fs. 

Besides  topics  such  as  bandwidth,  which  effect  each 
sampling  scheme,  there  are  other  concerns  that  need  to  be 
addressed  that  are  peculiar  to  the  particular  two-channel 
technique  chosen.  They  will  now  be  discussed. 

Conventional  Quadrature  Sampling 

Radar  and  communications  systems  are  becoming 
increasingly  dependent  upon  coherent  digital  processing. 
Conversion  of  signals  from  IF  analog  form  into  digital 
complex  samples  carrying  amplitude  and  phase  information  has 
been  t r ad i t i ona 1 1 y  implemented  in  the  form  of  two  parallel 
IF  to  baseband  converters  operated  in  quadrature  each 
followed  by  A/D  converters  which  thus  provide  digitized  in- 
phase,  p(t),  and  quadrature,  q(t),  components  (conventional 
quadrature  sampling  -  see  Figure  23).  Balancing  the  two 
baseband  converters  over  a  wide  dynamic  range  is  difficult 
and  phase  errors  are  typically  2°-3°  for  commercial  coherent 
detectors  (18:731).  To  examine  the  effects  of  phase  error 
in  the  quadrature  sampling  process,  recall  that  any  bandpass 


waveform  may  be  written 


x(t)  =  p(t)cosw0t  -  q(t)sinw0t  (42) 

Now,  Induce  a  phase  error,  0,  between  the  IF  to  baseband 
converters.  In  other  words,  assume  we  demodulate  with 
coswQt  and  sin(wot+0)  .  Beginning  with  the  cosine  channel 

x(t)coswQt  =  (p(t)cosw0t  -  q(t )slnw0t )cosw0t  (152) 

2 

x(t)coswQt  =  p(t)cos  w0t  -  q(t ) s inw0tcosw0t  (153) 

x(t)coswQt  =  ( 1/2 ) ( p ( t ) +p ( t ) cos 2w0t-q ( t ) s in2w0t ]  (154) 

The  low-pass  filter  removes  the  2w0  terms,  and 

x(t)cosw0t  ~  ( 1/2 ) p ( t )  (155) 

Similarly,  q(t)  is  obtained  from  the  sine  channel, 
x ( t )s in ( wot+0 )  =  p( t ) s in ( w0t+8 ) cosw0t 

-  q(t )sin(wot+0)sinwot  (156) 

But, 

sinAcosB  =  (1/2) (sin(A+B)  +  sin(A-B)]  (157) 

sinAsinB  =  ( 1/2 ) ( cos ( A-B )  -  cos(A+B))  (158) 

x ( t )sin ( wQt+0 )  =  ( 1/2 ) p( t ) ( s in ( 2wQt+0 )  +  sln0l 

-  ( 1/2 ) q ( t ) [ cos©  -  cos(2wot+0))  (159) 
The  low-pass  filter  removes  the  2w0  terms,  and 


x ( t ) s 1 n ( wq t +0 )  =  ( 1/2 ) [ p ( t ) 3 in6  -  q(t)cos9]  (160) 

Normally  we  expect  to  be  able  to  recover  the  amplitude, 
A(t),  and  phase,  #(t),  information  about  x(t)  using  the 


following  relations, 

A ( t )  e  [ p2 ( t )  +  q2 ( t  )  ] 1/2  (161) 

<*(t)  =  tan  ^  [  p  ( t  ) /q  ( t )  )  (162) 

What  we  actually  have  is 

A ( t )  ~  [p2(t)+(p(t)sin0)2+(q(t)cos9)2]1/2  (163) 

<*(t)  =  tan  1 { p ( t ) / [ p ( t ) s i n0  -  q(t)cos01)  (164) 


Obviously  the  phase  error  causes  attenuation  of  the  output 
signal.  For  small  fixed  errors,  this  is  tolerable.  For 
phase  errors  approaching  r90°,  however,  the  received  signal 
is  wiped  out  (17:206). 

A  further  cause  of  concern,  with  the  same  consequences, 
is  phase  errors  between  the  transmitter  and  the  receiver. 
Typically,  this  is  no  major  problem  when  the  transmitter  and 
receiver  are  in  close  proximity.  Again,  2°-3°  is  tolerable. 

The  key  point  to  thi3  discussion  is  that  random  phase 
errors  must  be  controlled  or  the  phase  errors  will  render 
this  conventional  quadrature  sampling  approach  useless. 

Quadrature  Sampling 

NOTE:  The  material  in  this  section  parallels  and  expands  on 


Persons’  quadrature  sampling  error  formula  (11:511-512). 

For  that  reason  and  ease  of  reference,  Persons'  notation 
will  be  used.  To  be  consistent  with  the  rest  of  this 
thesis  we  let  2nf0  =  w0  and  2nW  =  a. 

Recall  that  the  aim  of  quadrature  sampling  is  to 
recover  the  low-pass  components,  p(t)  and  q(t),  directly 
from  samples  of  both  the  bandpass  signal  x ( t )  and  its 
quarter -wavelength  translation  x ! t -1/ ( 4 f 0  )  ] ,  the  samples 
being  taken  at  the  low-pass  rate.  The  Grace-Pitt 
interpolation  formula  (7:1454). 

CD 

g ( t )  -  e  f (n/fs)s(t-n/fs) 
n  -  -  ® 

<T» 

♦  Z  f (n/fs+l/xf0)s(t-n/fs-l/xf05  (Iff: 
n  =  -  ® 

where  x  =  4,  allows  an  exact  reconstruction  of  high 
frequency  bandlimited  waveforms  with  a  minimum  average 
sampling  rate  of  W  samples/s  (or  W/2  samples/s /chan ne 1 ! 
whenever  2 fc  =  kw  as  in  Figure  29.  This  was  verified  1  y 
Persons  when  he  developed  a  quadrature  sampling  error 
formula . 

For  the  general  case,  where  2f0  does  not  equal  kw,  the 
interpolation  formula  given  by  Eq  (165)  with  x  -  4  will  not 

exactly  yield  f(t).  Persons  showed  a  way  to  measure  the 
amount  of  reconstruction  error  when  using  this  equation  for 
the  general  case.  More  spec  1 f 1 ca 1 1 y ,  he  determined  Q,  the 


Figure  29.  Bandpass  Signal  for  2£0 


kW  (k  =  integer) 


ratio  of  the  power  in  the  error  band  to  the  power  in  the 
signal  band. 

By  following  the  same  approach  taken  by  Persons,  it  13 
possible  to  more  fully  develop  the  quadrature  sampling  error 
formula  into  a  generalized  formula  in  which  the  effects 
on  Q  due  to  hardware  timing  errors  can  be  investigated. 

It  will  be  shown  that  timing  is  extremely  critical  when 
sampling  at  or  near  £s  =  o/n  samples/s  (average). 

Analy3 is .  Beginning  with  the  general  interpolation 
formula  of  Eq  (165),  recall  that  the  objective  is  to 
determine  the  power  in  the  error  relative  to  the  power  in 
£  (t)  . 

We  can  express  Eq  (165) 


g(t)  =  (f(t)E  o ( t -n/fg  )  ] *s ( t ) 


+  [f(t)E  <J(t-n/fs-l/xf0)  3*s(t) 


where  *  denotes  convolution. 

Taking  the  Fourier  transform  of  Eq  (166)  yields 


G(f)  =  fsS ( f ) { F ( f ) * IE  3(f-nfs)+E  fitf-nfg  )exp(-j2«nfs/x£0) ] } 


G  (  f  )  =  2fsS(f){F(f)  +  (1/2)E  F(f-nfs) [l*exp( - j«2nfg/xf0) ]) 

n^O 

(168 

Let  2 f s S ( f )  =1  in  the  region  where  F(£)  is  nonzero  and 
zero  elsewhere. 


2 f s S ( f  )  =  (recti  (f-f0)/W)  +  r ect [  ( f  +  f 0 ) / W ] ) 


This  i 3  equivalent  to  applying  an  ideal  antialiasing  filter 
to  the  bandpass  signal.  Solving  for  Sit) 


Sit)  -  ( l/2f„ ) (recti (f-fn)/W)  +  recti  ( f  +  fn)/wj )  (170 


The  Inverse  Fourier  transform  of  S ( f )  is 


s(t)  =  (W/fs)sinc(Wt )cos( 2«f0t ) 


Substituting  Eq  (170)  into  Eq  (168)  yields 


© 


G  (  £  ) 


=  F ( £ )  +  { (1/2) (recti (f-fQ)/w]  +  rect I ( £+£Q)/W] } 

L  F  (  n-nf  s  )  [l  +  exp(-zm2n  fs/xf0)  ]  )  (172 

n  /  0 

Notice  that  G(f)  is  made  up  of  F(f),  the  term  we  want,  plus 
seme  extra  terms.  It  is  the  power  in  these  extra  terms 
relative  to  that  in  F(f)  that  will  be  determined.  Because 
of  the  rect  functions  ( ant ia 1 i as  i  ng  filter),  only  the  terms 
giving  bands  in  the  vicinity  of  f0  will  interfere. 

Let  2fc/f£  =  z  +  r  =  K  ,  where  z  -  integer.  Als-, 
assume  that  f 3  >  W  ,  which  is  an  obvious  need  or  else  the 
repeated  spectral  bands  would  overlap.  Note  that  for 
2  £  -  '  f -  u  i  r.  t  e  g  e  r  ,  and  x  =  4  ,  there  is  only  c  n  e  band 

for  r.  =  i  c  which  overlaps  F!f),  but  the  amplitude  term 
1  +  exp: -cm)  -  C  .  In  general,  there  will  be  tw:  bands 
•over  lap  p  i  r.  g  Fif),  corresponding  t  •:  the  terns 
n  -  iz,  ♦  (  z+  1  )  .  Considering  only  the  posit; ve 

frequencies.  Figure  30  shows  the  z  and  !  z  ♦  1  )  term, 
overlapping  F  (  f  )  for  the  general  case  where  2£:  is  net  egu-. 
to  k.W  :  k  -  integer),  and  f  s  =  W  samples  s. 

Define  Q  equal  t .  the  rati  of  p  : we r  in  *  he  extra  band 
to  the  power  in  F ( f  1  ,  and  again  consider  only  positive 
f  r  e  g  u  e  n  c  i  e  s  . 

De  fine 

a  z  -  (1/2)11  *  exp!  j  n  4  f  c  z  >:  2  £  _  >  1  1  ”  3 

a  •  2  4  i  )  -  (  1  /  2  )  !  1  ♦  e  x  p  (  •  j  n  4  {  r  .  2  ‘  1  !  .  x  2  f  :  1  I  "  4 


*>t  si 


Figure  30.  Bandpass  Signal  for  2f0  j  kW  (k  =  integer) 


B-.  -  bandwidth  of  zth  term 


(175) 


E  r  2  ♦ 


bandwidth  of  (z+l)st  term 


(  176  ) 


Since  we've  applied  the  antialiasing  filter,  the  numbers  E 
ar'  *  ®  (  z  +  1  )  are  determined  by  the  overlap  wit!; 
rect  [  (  f-f0  )/w:  .  Let  the  power  in  F(f)  =  W  ;  then  from, 


JL 


«.n. 


.'tkWi 


Q  =  J  az I  ( Bz/W)  +  I  a ( Z  +  1 j  |  tB(z  +  1)/W] 


(177) 


where  Bz  and  B(Z  +  u  must  be  greater  than  zero,  else  they  are 
taken  as  zero. 

Now,  from  Eq  (173) 


|  a  z  |  2  =  ( 1/4 ) | £  1  +  exp(-j*4£sz/x2f0) ] j2 


(178) 


but  K  =  2f0/£s  and  z  =  IK]  ,  where  [  ]  is  the  greatest 

integer  function. 


| a z  j  2  =  ( 1/ 4  )  |  { 1  +  exp( —  j tc 4  IK  1/xK  ) }  |  2 


(179) 


|az | 2  =  {cos2 ( 2n [ K ] /xK ) } 


(180) 


Likewise, 


|  a  (  ^  + 1 )  |2  =  ( 1/4  )  |  { 1  +  exp{-jrt4(  (Kl+D/xK}}  |2 


(181) 


Is*  (  z  +  1 )  |  2  =  (cos2{2tc(  [ K ] +1 ) /xK }  } 


(182) 


Note  from  Figure  30,  the  center  of  the  zth  band  is  at 


£z  =  {zfs  -  f0}  =  { I K 1 2 f 0/K  -  f0) 


(183) 


and  the  center  of  the  (z+l)st  band  is  at 


£  ( z+ 1 )  =  Uz  +  l)fs  -  £01  =  {  ( ( K  )  +1 )  2£0/K  -  f0> 


(184) 


Now  determine  Bz/W  .  From  Figure  30,  Bz  =  fz  +  W  -  f0 


Bz/W  =  (fz  +  W  -  f0)/W 


Bz/W  =  £Z/W  +  1  -  £0/W  ( 

Bz/W  =  £Z/W  +  1/2  -  ( £0/W  -  1/2)  ( 

Substituting 

fz  =  [ K ] 2£0/K  -  f0  ( 

Bz/W  =  { [ K 1 2£0/K  -  f0)/W  +  1/2  -  {f0/W  -  1/2)  ( 

Since  f0  =  fsK/2  ;  K  =  2f0/fs 

Bz/W  =  1  -  { f s ( K- l K ) ) }/W  ( 

Likewise, 

B(z+1)/W  =  1  -  { f s ( 1-K+ [K I } /W  ( 


Substituting  into  Eq  (177)  and  using  the  definition 
K=z+r=  [ K 1  +r  and  L  =  £S/W 

Q ( K , L )  =  (cos2 ( 2n(K ]/xK ) } t [1-Lr  ]  ] 

+  {cos2 (2n(lK]+l)/xK)}{ [ 1-L( 1-r  )  ]  ]  ( 


where 


* 


I  Ll  U'U'U'U'U'Ll 


r  =  K  -  IK) 


L  =  fs/W 


(195) 


(196) 


and  l(  11  indicates  the  quantity  is  to  be  taken  as  zero  if 
it  is  negative.  Eq  (192)  holds  for  L  i  2f0/W  +1  ;  for 

larger  L,  the  error  is  zero. 

This  is  the  desired  result.  Note  that  substituting 
x  =  4  into  Eq  (192)  will  yield  the  quadrature  sampling 
error  formula  developed  by  Persons. 

Given  fQ  and  W,  Q  can  be  plotted  as  a  function  of  K  for 
any  L.  As  mentioned,  Eq  (192)  provides  us  with  the 
flexibility  to  investigate  the  effects  of  hardware  timing 
errors  on  Q.  If  we  use  Persons'  example  we  can  establish  it 
as  a  baseline  (he  assumed  no  timing  error).  In  other  words, 
we'll  assume  that  we've  sampled  £(t)  and  f(t-l/4f0)  with  no 


timing  error.  This  is  the  case  where  x  =  4 
L  =  f~/W  =1  we  have 


Q ( K , L+l )  =  {cos  (mz/2K)} ( [1-r] ] 


+  {cos  (u(z+l)/2K)}[[rl] 


Letting 


(197) 


which  is  plotted  in  Figure  31.  Figure  31a  shows  the  full 
scale  of  Q  for  0  <  Q  <  1.0  where  1.0  =  100%  error.  In  part 
b  the  Q  scale  was  expanded  to  emphasize  the  shape  of  the 
error  plot.  Notice  that  Q  is  zero  for  K  =  integer  and 
approaches  zero  for  narrowband  systems  (large  values  of  K). 


75 


The  interesting  question  is,  what  happens  to  Q  if  we 


have  a  timing  error?  In  other  words,  we  think  we’re 
sampling  f(t)  and  f(t-l/4f0),  but  we're  actually 
sampling  £(t)  and  f(t-l/xfQ)  where  x  does  not  equal  4. 
Since  quadrature  sampling  requires  sampling  the  bandpass 
signal  and  its  quarter-wavelength  translation,  we  expect 
some  reconstruction  error.  The  next  few  figures  will  show 
that  there  is  indeed  a  reconstruction  error  and  this  error 
can  be  quite  significant.  For  instance,  if  x=2  then  the 
bandpass  signal  and  its  half-wavelength  translation  are 
being  sampled  and  the  reconstruction  error  approaches  100% 
(the  signal  cannot  be  reconstructed). 

Q  is  plotted  in  Figure  32a  for  x  =  4.0,  3.9,  3.8,  3.7 
3.6,  3.5  and  0  <  Q  <  1.0  and  once  again  in  Figure  32b 
with  0  <  Q  s  0.16  .  The  individual  plots  for  the  various 
values  of  x  are  more  easily  distinguishable  in  part  b. 
Notice  as  the  timing  error  increases,  or  |x-4.0|  grows 
larger,  that  Q  increases.  Figure  33  shows  the  case  where 
the  timing  error  has  increased  to  the  point  where  we  are 
sampling  the  bandpass  waveform  and  its  half-wavelength  (x  = 
2.0).  In  this  figure,  Q  is  plotted  with  x  =  4.0,  3.8,  3.5 
3.0,  2.5,  2.0  .  As  before,  Q  approaches  100%  as  x 

approaches  2.0  .  Figure  33b  is  a  three-dimensional 

surface  plot  of  part  a. 

To  summarize,  consider  an  example  with  three  high 


frequency  bandpass  signals  with  center  frequencies  of  £2  - 
40  MHz,  £2  =  500  MHz,  and  £3=1  GHz.  Once  again 
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Q(K=2fo/W,fs=W=1)  vs  K 


TABLE  I 


Comparison  of  Timing  Errors  ( Tg  in  pS ) 


X 

Te  for 

fl  =  40  MHz 

f 2  =  500  MHz 

BBQ 

3.9 

160 

13 

6 

3.8 

329 

26 

13 

3.5 

893 

71 

35 

3.0 

2083 

167 

83 

2.5 

3750 

300 

150 

2.0 

6250 

500 

250 

fs  =  1.5W  ,  fs  =  2.0W  and  fs  =  5.0W  respectively.  As 

the  sampling  rate  grows  large  compared  to  the  bandwidth  the 
error  terms  approach  zero.  In  Figure  36,  Q  is  negligible 
for  x  =  4.0  (no  timing  error),  3.8  and  3.5  .  The  only 

significant  Q  terms  are  for  x  =  3.0,  2.5  and  2.0. 

Increasing  sampling  rates  will  reduce  the  Q  terms  even 
further . 

To  summarize,  quadrature  sampling  requires  very  precise 
timing.  Even  with  center  frequencies  less  than  100  MHz, 
nanosecond  precision  is  required.  If  such  precision  is  not 
possible,  higher  sampling  rates  will  be  needed  to  overcome 


the  reconstruction  error  terms. 


let  L  =  1  =  fs/W.  The  timing  errors,  TE,  is  a  function 
of  x : 

Te  =  T ( 1/4  -  1/x)  (198) 

where  T  is  the  period  of  the  corresponding  bandpass  signal. 

In  Figure  33,  a  value  of  x  =  2.0  results  in  a  Q  of 
almost  1.0  (or  100%)  .  Referring  to  Table  I,  and  putting 
this  in  terms  of  a  timing  error,  we  see  that  for  fj,  a  TE 
of  only  6nS  yields  a  Q  of  100%  .  Similarly,  TE  =  500  pS 

and  250  pS  for  f2  and  f3,  respectively. 

The  important  point  here  is  that  Q  is  directly 
proportional  to  the  center  frequency.  The  higher  the  center 
frequency,  the  more  precise  the  timing  must  be  to  keep  Q 
small . 

Typical  communication  signals  will  be  sampled  at  the  IF 
stages  of  the  receiver,  so  TE  for  f^  is  representative  of 
the  type  of  timing  precision  needed  for  quadrature  sampling 
when  fs  =  W  samples/s.  Referring  again  to  Figure  33  and 
looking  at  the  column  for  f]_  in  Table  I,  where  x  =  3.8, 

Te  =  329  pS  and  Q  is  relatively  small  for  large  values  of 
K.  A  value  of  x  =  3.0  represents  a  2  nanosecond  timing 
error  and  a  Q  that  levels  off  at  0.3  (30%).  Obviously, 
timing  must  be  kept  to  nanosecond  precision  or  it  won't  be 
possible  to  reconstruct  the  signal  from  the  samples. 

The  next  question  is,  what  can  be  done  to  get  around 
this  timing  problem?  The  answer  is  to  sample  at  a  higher 
rate.  The  next  three  figures  (34-36)  are  plots  of  0  for 


Q(K=2fo/W,L=fs/W=1 .5)  vs  K 


Figure  34. 


vs  K  with 


The  drawback  of  complex  sampling  Is  the  preprocessing 
that  must  be  done  to  obtain  the  Hilbert  transform  of  the 
bandpass  signal.  There  is  the  further  disadvantage  that  it 
is  not  easy  to  obtain  a  precise  Hilbert  transform  (the 
noncausality  of  the  Hilbert  transform  does  not  allow  one  to 
implement  it  for  signals  having  non-zero  bandwidth, 
theoretically).  Simpler  processing  is  usually  desirable 


(5:289)  . 


I v .  Conclusions  and  Recommendations 

Cone  1  us  Ions 

Five  bandpass  sampling  techniques  were  evaluated.  In 
each  case,  equations  for  the  bandpass  signal  were  developed 
in  terms  of  its  uniform  samples.  The  minimum  sampling  rate 
and  the  conditions  for  those  minimum  sampling  rates  were 
discussed.  In  comparing  sampling  rates  it  was  found  that 
each  of  the  two-channel  techniques  offered  a  rate  reduct  ior. 
of  up  to  one-half  over  the  s ingle-channel  approach.  There 
are  always  tradeoffs  that  result  in  varying  degrees  cf 
accuracy  and  must  be  considered  when  determining  wwi.;' 
technique  to  to  use.  For  example,  it  was  show:-,  tr.v 
conventional  quadrature  sampling  requires  pres..  ;■  ■ 
synchr on i za t i on,  quadrature  sampling  was  r-  »  * 

extremely  sensitive  to  timing  err::  ?  ar.  !  ; 

hard  to  implement  in  hardware.  It  w. 
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aspect  which  can  be  Investigated  is  the  reconstruction  error 
due  to  the  combined  effects  of  hardware  timing  errors  and 
varying  levels  of  noise. 
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The  purpose  of  this  thesis  is  to  investigate  and 
analyze  bandpass  sampling  techniques.  Bandpass  sampling 
theory  is  shown  to  be  an  extension  o£  the  Shannon  low-pass 
sampling  theory.  Five  bandpass  sampling  techniques  are 
investigated:  direct  sampling  (also  called  first-order 
sampling),  pth-order  sampling  (in  particular  second-order 
sampling  where  p=2),  conventional  quadrature  sampling, 
quadrature  sampling,  and  complex  sampling  (also  called 
Hilbert  transform  sampling). 

Direct  sampling  is  a  single  channel  approach  to 
bandpass  sampling  while  the  other  approaches  require  two 
channels.  The  two  channel  techniques  offer  a  sampling  rate 
reduction  of  up  to  one  half  over  direct  sampling,  but  there 
are  trade-offs  to  consider  when  choosing  a  two  channel 
technique.  It  is  shown  that  the  effects  of  random  phase 
error  can  render  the  conventional  quadrature  sampling 
approach  useless.  Timing  is  critical  with  quadrature 
sampling.  Reconstruction  error,  due  to  hardware  timing 
errors,  is  shorn  to  be  directly  proportional  to  the  center 
frequency  and  is  a  function  of  the  bandwidth  of  the  signal.  - 
Plots  were  developed  to  show  the  effect  of  timing  on  a  given 
bandpass  signal  when  the  center  frequency  and  the  bandwidth 
are  known. 
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